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Abstract 

We present an analytical model that enables throughput evaluation of Opportunistic Spectrum Or- 
thogonal Frequency Division Multiple Access (OS-OFDMA) networks. The core feature of the model, 
based on a discrete time Markov chain, is the consideration of different channel and subchannel allocation 
strategies under different Primary and Secondary user types, traffic and priority levels. The analytical 
model also assesses the impact of different spectrum sensing strategies on the throughput of OS-OFDMA 
J> ' network. The analysis applies to the IEEE 802.22 standard, to evaluate the impact of two-stage spectrum 

o ; 

sensing strategy and varying temporal activity of wireless microphones on the IEEE 802.22 throughput. 



Our study suggests that OS-OFDMA with subchannel notching and channel bonding could provide almost 



£ — ■ ten times higher throughput compared with the design without those options, when the activity and density 

of wireless microphones is very high. Furthermore, we confirm that OS-OFDMA implementation without 
subchannel notching, used in the IEEE 802.22, is able to support real-time and non-real-time quality of 
service classes, provided that wireless microphones temporal activity is moderate (with approximately 
one wireless microphone per 3,000 inhabitants with light urban population density and short duty cycles). 



Finally, two-stage spectrum sensing option improves OS-OFDMA throughput, provided that the length 
of spectrum sensing at every stage is optimized using our model. 
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I. Introduction 

One of the ways to combat artificial spectrum scarcity (2l is to augment existing radio access techniques 
with Opportunistic Spectrum Access Q (OSA). Wireless networks with OSA capabilities are able to 
search for unused licensed portions of the radio spectrum and communicate over those vacant radio 
frequencies whenever available radio capacity is insufficient, while meeting the required interference 
constraints. 

Orthogonal Frequency Division Multiple Access (OFDMA) is a multiple access technique where 
orthogonally-divided frequency subcarriers are assigned to individual users of the network. A subcarrier 
assignment is usually performed by a central entity, often Base Station (BS), and can be based on 
the quality of service (QoS) requirements of the individual users. Because of high spectral efficiency, 
as well as robustness against inter-symbol interference, OFDMA was a design choice of recent wireless 
networking standards, e.g. the IEEE 802.16 [4], the IEEE 802.20 © and 3GPP Long Term Evolution @|. 
As opportunistic spectrum use can be implemented efficiently with OFDMA, it seemed natural to connect 
the strengths of OFDMA with the flexibility of OSA. The first paper that introduced such concept, denoted 
in the reminder of this paper as Opportunistic Spectrum OFDMA (OS -OFDMA) was Q (referred therein 
as Spectrum Pooling) where OFDM subcarriers assigned to individual OSA users (denoted as Secondary 
Users (SUs)) are deactivated whenever the Primary User (PU) of the radio frequency band reappears. 
For a recent discussion on the topic of OS -OFDMA please refer to (H. 

So far no theoretical work on the system-level and cross-layer performance of OS-OFDMA networks 
has been reported. The need for theoretical framework for OS-OFDMA is also motivated by the recent 
introduction of OSA network standard IEEE 802.22 0-03] • The IEEE 802.22, an extension of the 
IEEE 802.16 standard, is designed to operate in the vacant TV band^j In the application domain the 
IEEE 802.22 has already been proposed to bridge remote wireless sensor networks with the command 
center lfl4l or support Internet connectivity in the rural areas 031 . 

Our goal is to develop the analytical framework to analyze the impact of traffic characteristics of 
OS-OFDMA network subscribers, the activity of the PUs of the radio spectrum, as well as the spectrum 
sensing algorithm and OFDM subcarrier assignment algorithm on the average throughput of OS-OFDMA 
network. Our approach is based on a cross-layer Markov chain-based analysis of OS-OFDMA, which 

'Note that the IEEE 802.22 is not the only networking standard proposed that focused on the operation in the TV white 
spaces. The remaining are recently published ECMA TC48-TG1 standard |12| focusing on porting local area networks to TV 
white spaces, and recently started IEEE 802.11af 1131 , similar in scope to the aforementioned ECMA activity. 



allows the investigation of interactions between medium access control and spectrum sensing layer. 
Since many options of OS-OFDMA subcarrier and subchannel assignment algorithms exist (namely, 
non-continuos subchannel assignment, as advocated by |7], and continuos version, as used in the IEEE 
802.22 J9l) it it is important to compare these designs using common analytical framework. In addition, 
analysis of two-stage spectrum sensing algorithm, proposed by the IEEE 802.22 standard, is challenging 
due to its complex effect on the medium access control layer and has not been explored in the context of 
OS-OFDMA communication. Finally, in the context of the IEEE 802.22 analysis could provide estimates 
of what QoS classes can be supported in OS-OFDMA, given realistic network conditions (such as number 
and type of primary users, the number and type of QoS classes enabled by OS-OFDMA network and 
the priority order in channel access for each class of users). 

Considering related work, in |fT6l a general framework of the IEEE 802.16 with OS A capabilities 
has been proposed with a very simplified networking model, based on Erlang-B formula |[T6l Sec. V- 
A], where the focus of the paper has been mostly on propagation calculations, including coverage, 
interference and protection distances. In [ 17 ] a simulation platform for the IEEE 802.22-like network, with 
limited set of ODFMA design options, has been presented. Focusing on the IEEE 802.22, interestingly, 
while many papers analyzed a certain functionality of the IEEE 802.22 network, like efficient spectrum 
sensing algorithm design ifTBl , lfT9l , circuit design for dedicated spectrum sensing EOl , MIMO extensions 
for the IEEE 802.22 physical layer Ell , game theoretic analysis of the the IEEE 802.22 networks 
coexistence 11221 - 11251 . 1126*1 . lf2"7l (with joint resource allocation), duplexing schemes lf28l (frequency 
hopping operation), |f"T31 (time division duplex design), and mesh establishment [29], it is desirable to 
develop unifying model that captures the IEEE 802.22 intrinsic features such as multiple classes of traffic, 
two stage spectrum sensing, different types of PUs and their temporal activity, and OFDMA subcarrier 
allocation process. 

The work closest to the scope of this paper can be found in lf30l in which the IEEE 802.16 system was 
evaluated. Obviously the model developed therein cannot be used directly to evaluate OS-OFDMA system 
due to the lack of spectrum sensing and PUs activity features. A work similar in our scope (analyzing the 
system level aspects of subchannel/subcarrier allocation strategies for OS-OFDMA) has been presented 
in |[3T1 . However, no comparison with the the IEEE 802.22 subchannel assignment has been considered. 
Furthermore, no QoS classes, PU priorities and two-stage spectrum sensing mechanisms were included 
in the model. 

Finally, we need to mention a set of papers analyzing performance of medium access control protocols 
for OS A networks. Some of the relevant ones include |32]-[35], however none of those works consider 
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Fig. 1. OS-OFDMA system model; SS: spectrum sensing. 

OFDMA, usually abstracting underlying physical channel structure. 

In this paper, we propose an analytical framework to quantitatively assess the performance of a net- 
work based on OS-OFDMA, considering features such as channelization structure, subcarrier allocation, 
resource assignment to network subscribers and different spectrum sensing methods. In the model we 
consider different priorities and channel dwell times of SUs and PUs of the spectrum. The developed 
model allows to calculate capacity (measured in terms of average throughput) for real-time and non-real 
time QoS classes of OS-OFDMA network. 

The rest of the paper is organized as follows. System model describing OS-OFDMA network design 
options in detail is presented in Section ITT] Analytical model for evaluating throughput of the considered 
system is presented in Section |III] Numerical results follow in Section [IV] Finally the paper is concluded 
in Section |\H 



II. System Model 

We consider a centrally controlled network with OS-OFDMA, where a BS manages resources and 
coordinates spectrum sensing of individual OS-OFDMA network subscribers. Although the proposed 
model is applicable to both uplink and downlink traffic, for simplicity we assume that only downlink 
traffic is transmitted. In this paper, we constrain ourselves to a single cell configuration with multiple 
SUs and multiple PUs, belonging to different user classes. This allows the exclusion of the effect co- 



channel and adjacent channel interference, as well as co-existence mechanisms in OSA network, on the 
investigation of the relation between PU type, its activity level and OS-OFDMA design options. On the 
other hand, we consider transmission errors due to fading and noise on the subchannels. 

A model of the OSA protocol stack under consideration is depicted in Fig. [TJ We identify two main 
modules: (i) data transmission, which is responsible for regular data communication between SUs, and 
(ii) spectrum sensing, which is responsible for efficient detection of spectrum opportunities for the OSA 
network; see also lf35ll for a similar model. Each component has its unique physical (PHY) and medium 
access control (MAC) layer. Obviously each layer has its unique design options, for example channel 
and subchannel management algorithms and multi-stage sensing. Also OSA network can be described 
by individual parameters such as type of traffic, activity level and bandwidth used. In Section III-AI we 
introduce specific OS-OFDMA system configurations considered in this paper. Later in Section Hl-BI we 
introduce the spectrum sensing and MAC design options of interest. We aim to calculate the average 
throughput obtained at the data MAC layer for all classes of SU traffic, which will be described in detail 
in Section UTAH 

A. OSA System Configuration 

1) Channel Setup and its Relation to OFDMA: The frequency domain consists of X channels, each of 
which is composed of Y OFDMA subchannels. The total number of subchannels is thus M = XY . Each 
subchannel is further composed of OFDM subcarriers. In this paper, we constrain ourselves to subchannel 
domain analysis. Furthermore, we assume that subchannel throughput C is on average constant, while 
its average value depends on the physical layer characteristics such as modulation, error control coding, 
and MAC layer overhead such as the OFDMA preamble length. 

In the time domain, transmission segments are divided in frames of length tt. At the beginning of each 
frame, SUs of OSA network detect the presence of the PU. We implicitly assume a synchrony between 
PU and SU activity, as it is a well established and classical assumption in the theoretical analysis, see for 
example recent publications of (33, (361 Fig. 2], (33 Sec. Ill], (371 Fig. 1], OH Fig. 2] ED, Il39ll-ll43l 
that follow the same path. Note, however, that the proposed model is extendable to the case where PU slots 
are offset in time from SU slots. This would require further analysis of PU channel access policies [44], 
|45l which is beyond the scope of this paper. We emphasize, as in |[35l , 11461 , that assumption about 
the synchrony between PU and SU connections allows to calculate throughput upper bounds compared 
to transmission on a slot-asynchronous interface. Furthermore, we assume that each node in the OSA 
network observes the same signal emitted by the PU, thus each SU performs all the sensing measurements 



individually and sends the measurement result to BS on the uplink. Then, BS makes a final decision 
about the presence of PU on each subchannel. 

2) PU and SU Types: We consider different types of PUs and SUs. For the PUs, depending on the 
bandwidth and the activity level, we classify them into: (i) a wideband PU (WPU) having low activity, 
i.e. with long busy and long idle times, and (ii) a narrowband PU (NPU) having high activity, i.e. with 
short busy and shorter than WPU idle times. This classification makes the analysis more detailed and 
realistic. It also makes different scenarios of interest possible to analyze. For example, WPU can represent 
wireless video links, while NPU can represent wireless microphones, both operating in the TV bands. 
For SUs, again making our framework general and applicable to multiple scenarios, we consider two 
types of users: (i) those receiving real-time traffic, denoted as the constant bit rate (CBR) SUs, and (ii) 
those receiving elastic traffic, denoted as variable bit rate (VBR) SUs, which are included in the the 
IEEE 802.22 standard. In our OSA network model, different types of SU traffic flows are generated at 
the upper layers, i.e. application, network and transport, and forwarded down to data PHY, while PU 
signals are detected at the spectrum sensing module, see Fig. [Q 

Furthermore, we assume that a hierarchical structure of users is present, such that the WPU has the 
highest priority in accessing bandwidth, NPU is the second in access hierarchy, followed by CBR and 
finally VBR. In other words, if users of different classes could access the same subchannel, the lower 
priority class user must vacate in order for the higher priority class user to utilize the subchannel. The 
evacuated CBR switches to the other idle subchannels or drops the connection if there is no idle subchannel 
available. For VBR, if the PU is detected, the active VBR connection squeezes the bandwidth [30, Sec. 
III-B] excluding the channel or subchannel occupied by the PU, and if there is no channel detected as 
idle, it buffers data until the PU disappears. Note that the behavior of VBR promises to obtain highest 
possible throughput, as demonstrated in |[35l Sec. V-B and Fig. 6], assuming no switching overhead, 
while CBR does not consider buffering due to the excessive delay that this class might experience while 
waiting for WPU or NPU to vacate the bandwidth. 

Because of the finite channel capacity, the number of users considered in the system is finite, but 
different for different user types. We assume that at most U w>max , f7 n ,max> U C:m3iX and U v>max of WPUs, 
NPUs, CBR and VBR connections, respectively, can be active at the same time in the considered 
bandwidth ||30l Sec. III-B]. For the data traffic of SUs and PUs, for analytical tractability, we assume that 
all users generate new connections according to the negative exponential distribution for the inter-arrival 
time and burst departure time, which is again a classical assumption in performance analysis studies |47l . 
The analysis can be extended to the general distributions of PU and SU traffic, which is beyond the scope 



of this paper. Note, however, that the recent measurement campaign [48 ] showed that more than 60% of 
measured PU activities distributions, including ISM and cellular bands, fitted an exponential distribution. 
The average inter-arrival and departure time are denoted, respectively, as 1/X W and l/fi w for WPU, 
1/A n and 1/ p, n for NPU, 1/A C and 1//U C for CBR, and 1/A^ and l/p, v for VBR. Also we assume 
that connection of each class except for VBR occupies a fixed number of subchannels. We denote the 
instantaneous number of subchannels utilized by a connection class as l w for WPU, l n for NPU, l c for 
CBR and l v for VBR. Note that the number of subchannels assigned to every connection is fixed and l c , 
l w , l n G N, except for a VBR connection. In that case l v € M., which stems from the fact that one data 
frame consists of a group of OFDMA symbols and the symbols in the OFDMA frame can be assigned to 
multiple VBR connections. Also, for VBR connections, the burst departure time depends on the number 
of subchannels used by VBR, thus 1//j, v is an average departure time when one subchannel is assigned 
to the VBR connection. 

B. Design Options 

1) Spectrum Sensing PHY and MAC Layers: Sensing PHY senses the PU signal and passes the 
measurement about subchannel availability to the sensing MAC layer for further processing. When a PU 
is present on the subchannel, it transmits a signal with a certain power and/or unique feature. Thus by 
detecting the power and/or the feature, the SU can decide whether the PU signal is on the subchannel or 
not. The main parameters for the spectrum sensing PHY are the probability of detection, the probability 
of false alarm and the sensing time. There is a trade-off between the sensing time and the resulting 
probabilities [49]. That is, if a SU takes a long time to sense a subchannel, the time for data transmission 
may be reduced. However, more idle subchannels can be detected because of high accuracy, which in 
turn may increase bandwidth utilization. Therefore, sensing time and sensing accuracy are the critical 
design options for the sensing PHY. 

In the sensing MAC layer, the SUs decide collectively, with the help of BS, on the PU state on the 
subchannel based on the sensing results of the sensing PHY layer. We denote the detection based on 
multiple users as collaborative sensing and that based on multiple periods as multi-stage sensing. Since 
the performance of collaborative sensing is relatively well known, see for example ll35l . in this paper, 
we focus mostly on multi-stage sensing} For the first results of multi-stage spectrum sensing in network 

2 We do not focus on recently proposed spectrum sensing methods for OFDMA networks, like |50| where quite-active sensing 
is proposed with non-active users sensing while others actively communicating, since they belong to a single-stage spectrum 
sensing category. 



context refer to 11511 -1531 for two-stage sensing multi-channel system, and [54] for two-stage sensing 
single channel system. 

In this work we limit ourselves to two-stage sensing, noting that our analysis can be directly extended 
to multi-stage sensing. The procedure works as follows. First, the SU senses the subchannel coarsely at 
the beginning of every frame, with short sensing time and low sensing accuracy. If the PU is detected, 
the SU switches to fine sensing mode (immediately, in the same frame), with long sensing time and high 
sensing accuracy. Depending on the sensing strategy, fine sensing can increase sensing accuracy ll54l or 
frequency resolution |5T| . Two-stage sensing can be described by different sensing PHY parameters for 
each stage. We denote the probability of detection as 5 a and 6f, the probability of false alarm as (fi a 
and cj)f, and the sensing time as r a and r/ for coarse and fine sensing stages, respectively. When setting 
8 a = 1, (f) a = 1 and r a = the two-stage sensing model reduces to a single stage sensing model. 

In order to evaluate the effect of spectrum sensing on the system throughput, we consider two unique 
sensing strategies. Firstly, we consider a sensing strategy where the SU senses all channels, including 
the operating channel, with coarse sensing and when the SU detects the PU on any of the channels, it 
immediately switches to fine sensing. We name this strategy as general two-stage sensing and denote it 
as So- 

Secondly, for a specific case when the bandwidth to transmit data is fixed and less than the whole 
allowed bandwidth, we investigate the following strategy. During coarse sensing the SU senses not the 
whole bandwidth but only a fixed bandwidth that is currently utilized for data transmission. If the PU 
is detected on the channel, the SU immediately senses all channels allowed to be utilized for the OSA 
system with fine sensing and switches to one of the the channels detected as idle. We name this strategy 
as two-stage active channel sensing and denote it as Si. Since in this strategy, in contrary to So, there is 
no need to always sense all channels, the sensing time is reduced. 

2) Data PHY and MAC Layer: Even though the OSA network is aware of subchannels being idle or 
busy, it should determine how to utilize the subchannels detected as idle for data transmission. There 
are numerous methods to utilize the idle subchannels in an OSA context, for a good overview we refer 
to |[35l . However, we selectively study four strategies that are proper for the centrally-controlled OS- 
OFDMA-based network. Those four strategies can be classified into two groups based on the purpose of 
the strategies. 

Firstly, we need to determine how much bandwidth is utilized for data transmission from the channels 
detected as idle. One strategy is to utilize all channels detected as idle from the allowed bandwidth. 
This strategy may maximize channel utilization at the cost of the wideband RF and signal processing. 



We name this strategy as variable channel/subchannel bonding and denote it as Bi. On the other hand, 
another strategy is to transmit data through only one channel, even though there may exist more channels 
detected as idle. Because the SU operates on the bandwidth of only one channel, the cost for the RF and 
signal processing is low. However, it is inefficient because some available bandwidth may not be utilized. 
Since it is one of the operating modes of the IEEE 802.22, also advocated by Federal Communications 
Commission |55], we also include it in our study. We name this strategy fixed channel selection and 
denote it symbolically as Bo- 

Secondly, we also need a strategy to avoid utilizing the subchannels on which the PU is detected. An 
efficient strategy is to notch out the subchannel detected as busy and utilize all other available subchannels 
for OS-OFDMA. We name this strategy subchannel notching and denote it as Ni. Note that we assume 
for simplicity that it is possible to notch out subchannels and transmit on the adjacent subchannels without 
causing interference, which is a common assumption in system level analysis, e.g. EH, BD1 . BTL [56 



On the opposite side, a conservative strategy is to exclude (block) all subchannels within the operating 
channel from accessing, even though only one subchannel is utilized by the PU, which is suggested for 
the IEEE 802.22. We name this strategy channel blocking and denote it symbolically as No- 

Please note that channel switching delays is in the order of tens to hundreds of microseconds ll60ll . lloTl 
which is a fraction of the channel sensing time. Dependent on protocol design, an additional OFDMA 
frame might be needed to communicate the decision about channel switch between base station and 
terminals. For simplicity channel switching delay is neglected in the analysis. 

3) Design Options of Interest: Because we have three groups of binary choices, i.e. S x , N x , and B x , 
where x € {0, 1}, there can be eight possible combinations of design options. However, not all options 
are feasible. First, we do not consider the combination of subchannel notching (Ni) and fixed channel 
selection (Bo), since it is a special case of N]Bi configuration with a single channel. Also, for two-stage 
active channel sensing (Si), we only consider channel blocking (No) and fixed channel selection (Bo), 
because two-stage active channel sensing (Si) is applicable to fixed bandwidth utilization case only. This 
leaves four combinations of the design options which are summarized in Table U Note, that all options 
except S0N1B1, are considered by the IEEE 802.22 standard. 

3 Please refer to a recent studies on that topic that prove the feasibility of such approach. For example, in 1571 a sidelobe 
suppression with guard band equal to only one OFDM subcarrier interval was shown. In |58| a OFDM subcarrier notching was 
proposed with only 4% of the available spectrum wastage. Finally a practical implementation of OFDM subcarrier suppression 
with perfect channel utilization at the cost of throughput reduction was demonstrated in 1591 . 
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TABLE I 

Summary of Design Options of Considered in The Paper 



Symbol 


Two-Stage Sensing 


Notching/Blocking 


Bonding/Separation 


S0N1B1 


General 


Subchannel Notching 


Subchannel Bonding 


SoNoBi 


General 


Channel Blocking 


Channel Bonding 


SoNoBo 


General 


Channel Blocking 


Fixed Channel Selection 


SiNqBq 


Active Channel 


Channel Blocking 


Fixed Channel Selection 



III. Proposed Analytical Model 

In this section we will describe the calculation of average throughput obtained using each of the 
considered OS-OFDMA designs. The analysis is based on a probabilistic framework utilizing Markov 
chain. We start with S0N1B1 option, which serves as a foundation to analyze the remaining three OS- 
OFDMA designs. 

A. Case SqN\Bi (General Sensing, Subchannel Notching, Subchannel Bonding) 

Generally, in OFDMA-based wireless networks throughput depends on how many subchannels are 
utilized in the idle spectrum by each SU connection type [30, Sec. III-B] (in case of our model: by 
every CBR and VBR connection). Let Y > T\(m c ,m v ,M a ) be the probability that m c and m v subchannels 
are utilized by CBR and VBR connections, respectiveljy when M a subchannels are detected as idle. 
In addition, let r)(M a ) be the average ratio of data transmission time to total frame length when M a 
subchannels are detected as idle. Then, the total system throughput, H, can be calculated as 

M 

H = Cj2v(M a ) Yl (m c + m v )Pn(m c ,m v ,M a ). (1) 

M a =0 m c ,m v e{0,- ,M) 

In addition, the throughput of CBR connection, H c , and VBR connection, H v , is computed using dT) 

by replacing m c + m v with m c for H c and m v for H v . In the subsequent sections we will describe a 

method to derive Pri(-) and r/(-). In this method we will hierarchically decompose Pri(-) and rj(-) into 

a set of conditional probabilities. Each probability will describe a particular relation between spectrum 

sensing outcome, state of CBR and VBR connection and the PU activity. 



In the paper we follow the convention that each newly introduced probability will be uniquely identified by a number and 
introduced with all argument variables, while its later callouts will be referred as Pr x (-). 
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1 ) Derivation of Probability of Number of CBR and VBR Connections and Available Subchannels, 
Pri(-); Values of m c and m v are easily determined if we know the number of CBR connections 
connected with the BS, U c , the number of VBR connections connected with the BS, U v , and the number 
of subchannels detected as idle, M a . Thus 

t/c.max C/„,m„x 

Pr i(')= £ J2 Pr 2 (m c ,m v \U c ,U v ,M a )Pr 3 (U c ,U v ,M a ), (2) 

u c =o u v =o 

where Pr2(-) denotes a set of allowed subchannel configurations occupied by CBR and VBR connections 
and Prs(-) denotes the probability of active U c CBR connections, U v VBR connections and M a idle 
subchannels. 

a) Derivation of Allowed CBR and VBR Subchannel Configurations, Pr2(-); Since the total number 
of subchannels used by all CBR users, UJ C , cannot be greater than the number of the available subchan- 
nels, M a , because if there is no subchannel available the CBR connection will be blocked, the only valid 
case is U c l c < M a . Furthermore, a CBR connection has a higher priority than a VBR connection, all 
CBR connections can transmit data through all m c = U c l c subchannels. Then the remaining subchannels, 
i.e. m v = M a — U c l c , are used by VBR connections. On the other hand, if there are no VBR connections 
in the system then m v = 0. Therefore, defining U(x), where U(x < 0) = and U(x > 0) = 1 we have 



Pr 2 (-) 



1, U c l c < M a , m c = UJ C , m v = (M a - UJ C )U(U V ); 

w ) 

0, otherwise- 



b) Derivation of Probability of Active CBR and VBR Connections and Idle Subchannels, Pr3(-): 
Since our model considers arrival process of SU connections, which departure process is affected by PU 
temporal activity, to calculate Prs(-) in (O, we need a tool to evaluate steady state probability of given 
number CBR and VBR connections, as well as the number of idle subchannels. To do this we use a 
widely used method based on the composition of a Markov chain |47l Ch. 11]. 

We introduce a Markov chain state {U c , U v , M a }. Furthermore, we introduce the state transition 
probability which describes the change in {U c ,U v ,M a } between time t — 1 and t, denoted as Pr4 (•). 
Then we can compute Pr3(-) by solving the Markov chain, given X]x Pr3 (') = ^ anc ^ Pr 3(') = 
SjVf(t-i) P r 3(') P r 4(')> where A4 is a set of all possible states {U c ,U v ,M a } and My°> is a set of the 
states at time x. Based on the conditional probability property and independency of the variables, we 
decompose Pr4(-) as 

Pr 4 (.) ^ Pr 5 [ui*\U<*)\ut 1 \ut 1) ,M^,Mt- 1) ) ^ (m^M^) , (4) 
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where Prs(-), denotes the probability of change in number of CBR and VBR connections count, and 
Ptq(-), denotes the probability of change in number of subchannels detected as idle. We derive those 
expressions below. 

c) Derivation of Probability of Change in CBR and VBR Connection Count, Prs(-): We can 
decompose Pr5(-) into probabilities denoting a change in the number of connections separately for 
CBR and VBR. Accordingly, we define Pr 7 (u^ \U^~ l \ M%\ M^~ 1} ) as the conditional probability 



of the number of the CBRs at time t for the given number of subchannels detected as idle, and 
Pr 8 ( UP |[4* _1) , UP , C/ c ( * _1) , M® , M^ 1] ) as the conditional probability of the number of VBR con- 



nections at time t for the given number of CBR connections and the available subchannels. Note that 
in Pr7(-), since we assume that the CBR has higher priority than the VBR and the VBR connections 
utilize the remaining subchannels after subchannel assignment for all CBR connections, there is no 
dependency on [/] . Furthermore, note that in Pr7(-) and Prg(-) change in the number of CBR and 
VBR connections, respectively, depends on the number of subchannels detected as idle at time £ — 1 and 
t. Then we have 

Pr 5 (•) ^ Pr 7 (uMlUt^MPM*-*) ^8 (u^\Ut l) ^\U^ l \M^ ,M^ l) ) ■ (5) 

We proceed with describing the process of deriving the expressions for Pr7(-) and Prg(-). 

d) Derivation of Probability of Change in CBR Connection Count, Pry(-): First, we consider valid 
conditions for U c , Uc , M a , and M a for Pr7(-) involving all possible number of i arriving and 
j departing connections. We denote the number of users being able to utilize all available subchannels as 
U { a ] = \M^/l c \ at time t and U^ 1] = [M^ 1] /l c \ at time t - 1. Then, because U^ , U^~ 1] denote 
the maximum number of users Uc < U a and Uc < Us, • Furthermore, because the possible sets 
of i and j are different for the cases Uc < Ua, U c = U a > 0, and Uc = U a = 0, we consider 
them separately. 

The first case, U c < U a , represents the situation when the number of subchannels detected as idle 
is more than the number of all subchannels that will be utilized by CBR connections before spectrum 
sensing. In other words, no CBR connection is blocked due to the PU appearance. Because the number 
of CBR connections is Uc at time t — 1 and U c at time t, the change in the number of CBR 
connections is i — j = Uc — Uc ■ In addition, because there are Uc active connections at time 
t — 1, more than Uc connections cannot be released, i.e. j < Uc ■ Therefore {i,j} € fC c , a — 

{i,j\i-j = u^-ut 1 \j<ut 1) }. 

The second case, Uc = U a > 0, denotes the situation when CBR connections may be blocked 
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due to the PU arrival. Then, before spectrum sensing, the total number of connections including newly 
generated connections is Uc + i — j. However, after spectrum sensing, the connections that utilize 
subchannels detected as busy are blocked and the remaining connections, Uc , utilize all available 

h 



subchannels, Ma . Thus, the number of CBR connections before spectrum sensing, Uc +i 



can be greater than or equal to the number of CBR connections after spectrum sensing, Uc , but should 



be less than or equal to U cmax , i.e. Uc < Uc + i — j < U cmax . 



Therefore {i,j} € K. 



c.b 



{'■ 



m 



w 



U? 1} <i-j< Uc. 



rW 



At) 



The third and final case is when Ua = Uc = 0. In this situation, because there is no subchannels 
available for CBR, the number of CBR connections should also be zero. Consequently, regardless of i 
and j, the conditional probability Pr7(-) under this condition is always one. 

Now, we introduce two supporting functions, G x (i\U x , X x ) and T x (j\U x , jx x ), which will be used to 
derive Pr7(-), where G x (i\U x , X x ) is the probability that i connections are newly generated from U x 
available users with arrival rate X x , and T x (j\U x , fj, x ) is the probability that j connections are released, 
each with departure rate \i x . G x (i\U x , X x ) and T x (j\U x ,/i x ) are derived in Appendix lAl Our approach 
to derive Pr7(-) is to calculate all possible sets for i and j and applying them to G x (i\U x , X x ) and 
Tx(j\U x ,/j, x ). In result, Pf7(-) is derived as 



Pr 7 (-) 



(*) 



Z {l , 3} eic^ T cUPr^c)Gc(i\U 



r(*-l) 



A, 



rr(*-i) < rr(*-i) 
U® < U^; 



At) 



Ewjeic, TctiPc . ^)G C U\ur L > +i-j,x 



At- 1 ) 



(6) 



1, 



r/(*-i) < rr(*-i) 

UP = U® > 0; 
rr(*-i) < ir(t-i) 

crw = u® = o. 

e) Derivation of Probability of Change in VBR Connection Count, Prs(-); In our model we do 
not consider the case that the VBR connection is blocked by the PU because VBR connections are 
assumed to be buffered instead of blocked when there is no available subchannel. Thus, assuming that all 
VBR connections share the same portion of the idle bandwidth, we calculate the number of subchannels 
assigned to one VBR connection, l v , as 



'v 



Af(t-i)_ c/c*- 1 ) 



Ul, 



(t-1) 



ut 1] > 0, 



0, 



U 



(t-1) 



(7) 



0. 
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In turn, Pr§(-) from (f5]) is denned similarly to ((6]) as 

Pr 8 (-)= Yl TvVluPXltvlGv^luV-V,)^), (8) 

{i,i}e/c„ 

where K v = U,j\i-j = uP - U^ 1] ,j < U^' 1] \. 

f) Derivation of Probability of Change of Subchannel Count Detected as Idle, Prg(-). - Proceeding 
to derive Prg(-) in (01), it can be decomposed as 

Pr 9 (mPM-'A 

Pr 6 (-) = V 7 , t 1)N J , (9) 

Prio [Mi X) 



where Prg(-) denotes the probability that M a subchannels were detected as idle at time t — 1 and t, and 
Prio(-) denotes the probability that M a subchannels were detected as idle at time t — 1. We will explain 
the derivation of Prg(-) and Prio(-) in subsequent sections. 

g) Derivation of Probability of M a Number of Subchannels Detected as Idle at Time t—1 and 
t, Prg(-); The idea behind derivation of ( fTOl ) is that number of detected subchannels depends on what 
sensing stage was utilized at time slots t — 1 and t and how many NPUs and WPUs were present at these 
time slots. It can be defined as 

pr 9 (.) ^ Yi p ni (^^r^ur^^^v^- 1 )) 

j7(t-i) {/(t-i) £(t-i) 

xPr 12 ([/w 1 e i) 1 ^ , ,e i u (<) ,^) ) do) 

where S^ and S^*" 1 ) are the sensing stages at times t and t — 1, respectively. Furthermore, Prn (•) 
denotes the probability of the number of subchannels detected as idle given the number of WPUs and 
NPUs and the sensing stage at time t-1 and t, and Pri 2 (u$ , U^~ 1] , u£\ U%~ 1] , S^\ S^A denotes 
the joint probability of the numbers of WPUs and NPUs and the sensing stage at time t — 1 and t. We 
can further decompose Prn(-) and Pri2(-) into subsequent probabilities. 

h) Derivation of Probability of Number of Subchannels Detected as Idle Given the Number of 
WPUs and NPUs and the Sensing Stage at Time t — 1 and t, Prn(-); Probability Prn(-) in (fTOl) can be 
decomposed into products of probabilities describing available number of subchannels detected as idle 
at time t and time t — 1, because the number of subchannels detected as idle and the state of the sensing 
stage at time slot t is independent from time t — 1. That is 

Prn (•) 4 Pr 13 (il^*- 1 ^*- 1 ),^ 1 ),^ 1 )) Pr 13 (mM\U®,U®,sW) , (11) 

where Pri3(-) is derived in Appendix iBl 
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i) Derivation of Probability of the Numbers of WPUs and NPUs and the Sensing Stage at Time t — 1 
and t, Pri2(-).' Probability Pri2 (•) in (fTOb can be decomposed into conditional probabilities as follows 

PnaO ^ Pr 15 (S®\U$,U®) Pr 15 (s^^U^, U^) 

x Pr 16 (u$\ljV\Ut l \ut l) ) Pnr (u^M^) , (12) 

where Pris(-) is the probability of being in a certain sensing stage given the number of NPUs and WPUs, 
which is derived in Appendix ICl and Prig(-) and Pri7(-) are the state transition probability and the steady 
state probability for the state {U w ,U n }, respectively. As in the case of Prs(-) denoting probability of 
active number of CBR and VBR connections, since the state of NPU and WPU changes randomly and 
independently from time slot to time slot, we can apply the same method of Markov chain construction 
to derive the probability of being in any of the {U w , U n } states. We will describe their derivation below. 
j) Derivation of Steady State Transition Probability for State {U w , U n }, Pri7(-): Probability Pri7(-) 
in (fT2l ) is computed by solving a Markov chain, such that J2u P r i7(') = 1 an ^ P p i7(') = X^c-d P r i7(') P r i6( - )> 
where U is the set of all possible values of U w and U n , and U^ 1 ^ is the set of the same parameters 
at time t — 1. Based on the assumption that the WPU has a higher priority of channel access than the 
NPU, the state transition probability Prig(-) is derived as 

Pr 16 (•) 4 Pr 18 (u^\Ut 1] ) Pri 9 (u^U^, U$\ U^) , (13) 

where Pris(-) denotes the probability of change in number of WPUs, U w , between time slot t — 1 and 
t, and Prig(-) denotes the probability of change in number of NPUs, U n , between time slot t — 1 and t. 
We describe their derivation below. 

k) Derivation of Probability of Change in the Number of WPUs and NPUs, Pris(-) and Prig(-): 
From the introduced model, the inter-arrival time and departure time follows the negative exponential 
distribution. Therefore we can use equations (l24l ) and d25l ), derived in Appendix[A] in the similar way as in 
the derivation of © and ([8]). By denoting the available subchannels for NPU as U e 2 = (M — Uw l w )/l n 



at time t and U, 



el 



(M — Uw lv))/l n I at time t — 1, we can derive Pris(-) as 



Pr 18 (-)^ J2 T w (j\Ut 1] ,^)G w (i\U^,xJj, (14) 

{i,j}e/c w 

and Prig(-) as 



Prig(-) 



J2 T n {j\u!£\ii n )G n (i\ut 1] An) , Ut l) < U el ,U® < U e2 ; 

{i,j}eK n , a 

J2 T n {j\U^^ n )G n {i\Ut 1] + i~JAn), Ut X) <U el ,uV = U e2 , 



(15) 
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where JC W = U,j\i - j = U$ - U^~ X \j < % |, /C„ >a = U,j\i - j = £/„ - f/A , 
j < Ut 1] } and K nfi = {i, j|C/W - ET^" 1 ) < j - j < [/ njmax - C/A* _1) , j < C# -1) }. 

Zj Derivation of Probability of M a Subchannels were Detected as Idle at Time t — 1, Prio(-)-' 
Finally, Prio(-) in (0 is calculated as 

Vt/i* -1 ',!/, 1 ;'" 1 ',^'*- 1 ) 

x Pns^- 1 )!^- 1 ),^- 1 ))?^^- 1 ),^- 1 )). (16) 

2 J Derivation of Sensing Overhead, r/(-): The value of ry(-) depends on the sensing stage since a 
longer sensing time for one stage can reduce rj(-), while a shorter sensing time for another stage can 
increase ■>](■). We recall a variable S, introduced in Appendix |B1 and Appendix ICl which indicates the 
sensing stage, such that S = denotes the case when the OSA network performs only coarse sensing 
without switching to fine sensing, and S > denotes the case when the OSA network performs coarse 
sensing and switches to fine sensing immediately. Specifically, 5 = 1 denotes the case when the OSA 
network detects the idle subchannel and 5 = 2 denotes the case that no idle subchannel is detected so 
that the network waits until the next sensing period without transmitting data. 

Defining Pro (5, M a ) as the joint probability that the current sensing stage equals to 5 and the number 
of subchannels detected as idle is M a , we can compute i] as 

r,(M a ) = Pr o (0 ; M a )^f^+Pr (l, M a ) tf ~ T * +Tf . (17) 

Note that there is no Pro (2, M a ) in (TTTT ) because for 5 = 2 no data can be transmitted, so the ratio of 
the data transmission time to the total frame length is zero. Finally, we can derive Pro(-) in (fTTT ) as (fT6l ) 
removing S^^ 1 ' from the lowest bound of summation. 

B. Case SqNqBi (General Sensing, Channel Blocking, Channel Bonding) 

There are two major changes for the analysis of SoNoBi case in comparison to S0N1B1. First, the 
number of subchannels detected as idle should be the integer multiples of Y because in this case the 
smallest quantity of idle bandwidth is one channel. Second, the number of available subchannels is 
determined not only by the number of the NPUs but also by the position of the NPUs in the spectrum. 
For example, if two NPUs appear on different subchannels in the same channel, the SU cannot utilize 
that channel. If NPUs occupy subchannels located in two different channels, those two channels cannot 
be used by SUs. Considering those changes, we need to modify probabilities related to the number of 
subchannels detected as idle and the sensing stage used, that is Prn(-) in (fTTT) . Pri3(-) in ( |28T ) and Pris(-) 
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in d29l ). Note that the other probabilities remain the same as derived in Section IIII-Ai In the process of 
modification of the above expressions we assume 5f ~ 5 a ~ 1 to reduce the complexity of calculations. 
In general, the OS-OFDMA system needs to keep a high detection probability to protect the PUs, which 
makes this approximation reasonable. Note, however, that we will still consider the effect of false alarms. 
Thus, even if there is no PU on the spectrum, the SU may falsely detect an idle subchannel as busy. 
Moreover, without loss of generality, we assume that a WPU occupies one channel, i.e., l w = Y. 

To modify Prn(-) from (TTTb we observe that the change of the number of PUs between time t and 
t — 1 may affect the change of the position of the PUs and, as a result, can make an impact on the number 
of channels detected as idle at time t. Thus, introducing the conditional probability of the number of 
idle subchannels as Pr 23 (m^\U^\u^\ S®, M { a~ l \ U^~ 1] , U^~ l) , S^A, Pr n (-) in CD} is newly 
defined as 

Pr 11 (-)^Pr 13 (Mi*- 1 )|t4f- 1 ),^- 1 ),5( 4 - 1 )) 

x Pr 23 (Mi\Mt l \U^\ut l \u!t\ut l \S^\S^) , (18) 

where Pr23(-) is derived in Appendix [E] Finally, Pri3(«) and Pris(-) are modified in Appendix |Pl and 
Appendix |F] respectively. 

C. Case SqNqBq (General Sensing, Channel Blocking, Fixed Channel Selection) 

The analysis in this section is based on the derivation from Section ITlI-B I since this option still considers 
subchannel non-notching. The major change here is that we perform analysis for data transmission on 
one channel instead of all X channels. This change affects the valid condition for Pr7(-) in Q and 
l v in (O. Because the maximum number of available subchannels is limited to Y, Ua and U a 



becomes 



r(*) 



r(t-i) 



mm[Y/l n ,Ma /l n ) and min I Y/l n ,Ma /In) , respectively. Moreover, considering that 



the maximum number of channels utilized by VBR is also limited to one channel, then (0 should be 
modified as 

W* ' ^ >0 < (19) 

o, ut 1] = o. 

Also, we need to modify Prg(-) in © (probability of change of subchannel count detected as idle) 
considering the limitations of the available subchannels. Even if a SU detects more than one idle channel, 
the SU will utilize only one channel. In terms of the definition, we have to sum all probabilities that the 
number of subchannels detected as idle is greater than or equal to Y, in order to compute the probability 
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that one channel, i.e. Y subchannels, are detected as idle. Thus, 

Pr 9 (0,0) 
Prio(O) ' 



Pr 6 (- 



Ef =y Pr a (0,s) 

Ef =y Pr 9 Qr,0) 
Prio(O) ' 

H*. yeiY ,.., M} rr B (x,y) M (t) > Q M (t-1) > Q 



Mi*) 


= 


M ( a~ 


-1) 


= 


M^ 


= 


M { t 


-1) 


>0 


m« 


>0 


Mt 


-1) 


= 



(20) 



Ef^PnoOr) 

D. Case S\NqBq (Active Channel Sensing, Channel Blocking, Fixed Channel Selection) 

In this case the SU performs the coarse sensing for only one channel currently utilized for data 
transmission. Thus, the PU on a channel not used by the SU does not affect the OSA network, and as 
a result, the sensing stage of the OSA network becomes more sensitive to the location of the PUs in 
the radio spectrum. Therefore, probabilities related to the sensing stage S, such as Pri2(-) of (fT2l and 
P r i5(') of (1291 in Appendix ICl and Pr2o(-) of (fJTJ) in Appendix |Pl need to be updated. 

First, we update the definition of Pr2o(-)- m option SiNoBo, even if more than one channel is detected 
as idle, only one channel is utilized for data transmission. Thus, we consider only the case when one 
channel was detected as idle even though there can be more channels detected as idle. Then, obviously 
if a OSA network detects an idle channel, i.e. S < 2, the number of channels detected as idle is one. If 
S = 2, the probability that no idle channel is detected is also one. Thus Pr2o( - ) can be modified as 

a (l, X a = 0,S = 2orX a = l,S<2, 
Pr 20 (-) = < (21) 

I 0, otherwise. 

Next, we present the modification of Pri2(-) in (fT2l . This modification is based on the fact that the 
sensing stage at time t can be affected by the number and the position of WPUs and the NPUs at time 
t - 1. Thus, denoting Pr 25 (^S®\U$,UJP,S ( - t ~ 1 \u£~ 1 \u£~ 1) ') as the conditional probability of the 
sensing stage at time t given the number of NPUs and WPUs for times t and t — 1, Pri2(-) is modified 
as follows 

Pr 12 (-) ^ Pr 15 (sP-Vpt^Ut**) Pne (v^^Pt^Vt^) 

x Pr 17 (ut'Kut^) Pr 25 (S®p£\ E#\ S^M^ M^) • (22) 



For the same reason as in case of Pri2(-), Pris(-) in (1291 ) needs to be modified as well. The modification 
process is presented in Appendix [Gj 
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IV. Numerical Results 

In this section we provide performance results for all considered options. We note that all analytical 
results were verified via simulations using a method of batch means for 90% confidence interval. Each 
simulation run, with a warm-up period of 10000 network events, was divided into 100 batches, where 
each batch contained 10000 network events. 

Due to many parameters considered we limit our numerical investigation to three most representative 
case studies. That is, we consider the impact of varying number of NPUs, the impact of varying PU 
activity and the impact on two-stage spectrum sensing design on the system throughput. Results are 
presented in Section IIV-BI Section IIV-CI and Section IIV-D[ respectively. 

A. Calculation of Average Subchannel Capacity C 

Before we proceed with the presentation we need to comment on the calculation of C. Following the 
IEEE 802.22 model @ we consider the PHY capacity C p and MAC layer overhead £ separately, such 
that C = (1 - C)C P , where C p = 460.8 kbpsfi 

For the MAC layer overhead £ we consider the frame structure of the IEEE 802.22 such that one 
downlink OFDM symbol of all subchannels is assigned to a preamble and two downlink OFDMA symbols 
of all subchannels are assigned to management messages. Also considering errors on the subchannels, 
we assume a bit error rate of 10~ 6 . Hence, we calculate the MAC layer overhead reduction factor as 
1 — £ = 0.8125. Therefore the total subchannel capacity is C = 374.4 kbps. 

B. Impact of Varying Number of NPUs on OS-OFDMA Design 

Throughout this section we assume that X = 4 channels are available and split into Y = 10 
subchannels. Frame length is set to tf = 20 ms, which represents the length of two frames of the 
IEEE 802.22 [9], i.e. it represents the time SU device searches for OFDM preambles on a given channel 
and highly conservative value of a inter-frame sensing interval j9j Table 233]. False alarm and detection 
probability for the coarse sensing case is 5 a = 0.99 and <j) a = 0.1, respectively, while for the fine sensing 
case 5 a = 1 and <fi a = 0, respectively. Sensing time during coarse sensing is r a = 0. Since we have 
assumed that uplink and downlink are divided by time division duplex and coarse sensing is performed 

5 Assumptions: 16-QAM modulation with 4 bits per OFDM symbol and 1/2 channel coding per subcarrier; for guard band 
for NPUs, six subchannels of the IEEE 802.22 represent one subchannel in the system model. Uplink and downlink are time 
division duplexed, where 16 symbols and 8 symbols are assigned to downlink and uplink, respectively. 
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in the uplink, as in O, the sensing overhead is zero. In the fine sensing phase tj = |iy, modeling 3 
consecutive frames used for the fine sensing and the following 2 frames used for data transmission until 
next sensing period. Furthermore we assume that l c = 1, l w = 10, l n = 1 subchannels are allocated to 
CBR, NPU and WPU connection, respectively. Note that the number of subchannels assigned to VBR 
connection, l v , is varying depending on the network state. The maximum number of users of each class is 
U v>msx . = 2, U CjmsK = 10 and U w , max = 2 (investigating impact of NPUs), and U v>max = 2, £7 c , max = 10 
and U rhmSiX = 10 (investigating impact of WPUs). Those values correspond to a small network and allow 
for an easier understanding of the subsequent numerical results. 

Also, for the purpose of this section we set up inter-arrival and departure times taking into consideration 
the IEEE 802.22 network, where many active licensed users operate over the TV band. Since in general, 
the traffic pattern of PUs for such case is not well known (more discussion on this aspect is presented 
in Section IIV-CI) . for the WPU we keep the wireless assist video devices J62l in mind, which can be 
assumed to broadcast on average four hours of signal transmission for every twelve hours on average in 
this scenario. For the NPU, we consider environment with numerous wireless microphones and assume 
that they appear every two hours on average and utilize channels for one hour on average. For CBR, 
considering voice or video transmission, we assume that on average five minute long CBR connection 
is generated for every five minutes on average. For VBR, we assume that data traffic is generated every 
two hours on average and continues for two hours if one channel is assigned for a VBR connection. 
For simulation efficiency, since we operate in large parameter ranges, we scale them down by setting 
the CBR connection arrival rate to one second with preserving the ratios between all traffic parameters, 
i.e. we normalize average inter-arrival and departure times of all users in the unit of five minutes by 
dividing them by 300 seconds. Thus for the large number of users, the inter-arrival time becomes shorter. 
For the analysis, we calculate the inter-arrival time by dividing the individual inter-arrival time by the 
maximum number of users. In summary, l/X w = 144/{7 u , jmax s, 1/A„ = 24/C/ nimax s, 1/A C = l/£/ c>max 
s, 1/Xy = 12/U v>max s, l/i2 w = 48 s, l/fi n = 12 s, l/// c = 1 s and ljfi v = 240 s. 

The results are presented in Fig. |2] The throughput of every OS-OFDMA design option decreases 
with increasing numbers of NPUs and WPUs. First, we observe that S0N1B1 is the best design option 
when total and VBR throughput is concerned, which is due to the highest flexibility in exploiting all 
spectrum opportunities. Second, interestingly with low number of NPUs and WPUs CBR throughput is 
higher for SiNoBo than for S0N1B1, while with the high number of WPU and NPU the opposite holds, 
see Fig. |2(b)| and Fig. 2(e) This is because of the sensing overhead of the fine sensing that is performed 



more frequently for SqNiBi than for SiNqBo- Third, there is no difference in CBR throughput between 
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Fig. 2. Performance of different OS-OFDMA options as a function of NPU, Fig. (a), (b), (c), and WPU, Fig. (d), (e), (f) 
for X = 4, Y = 10, M = 40, t s = 20ms, 5 a = 0.99, <j> a = 0.1, r a = 0, 6f = 1, <f>f = 0, r s = ft/, C = 374.4kbps, 



(7 m ,max = 2 (WPU), (7 n , max = 10 (NPU), L = 10, « n = 1, C/ C , max = 10, l c = 1, (7„, ; 
(b) (e) throughput of CBR connections, and (c) (f) throughput of VBR connections. 



= 2; (a) (d) total system throughput, 



SoNrjBi and SrjNoBo when the number of NPUs varies. This is because in this network setup CBR, which 
has higher priority than the VBR, utilizes enough resources even though channel bonding is not applied. 
Fourth, the SoNoBi design option is extremely sensitive to the activity of the NPUs. The throughput of 
this design option decreases rapidly as the number of NPUs increases, compare Fig. |2(a)| with Fig. |2(c)| 
This effect is not visible however when the NPU number is kept fixed, but the number of WPUs changes, 
see Fig. |2(d)| and Fig. 2(f) and compare with Fig. 2(a) and Fig. |2(c)[ respectively. This is because SoNoBi 
is sensitive to the position of NPUs in the spectrum, i.e. the lack of subchannel notching (No), causes 
this option to perform worse. 

Further, we demonstrate a benefit of active channel sensing strategy in two-stage sensing. Surprisingly, 
the total throughput of S1N0B0 is greater than that of S0N0B1 for large number of NPUs and WPUs, see 
Fig. 2(a) and Fig. |2(d)[ even though the implementation for S0N0B1 enables wider bandwidth sensing 
than SiNqBo- This is because for SoNqBi, due to channel bonding, the probability that OSA network 
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needs to sense the channel is much higher and network needs to perform sensing often. While for S1N0B0 
only actively used channel is sensed. 

Note that the simulation results agree with the analysis in all figures. A slight mismatch between 



simulation and analysis for cases S0N0B1, SoNoBo and S1N0B0 in Fig. |2(e)| is a result of the approximation 
used in calculating throughput for these design options, see (l36l ). (l40l ). More discussion on this issue is 
presented in Section IIV-CI 

C. Impact of PU Activity on OS-OFDMA Design 

To investigate the impact of varying PU activity on the performance of different OS-OFDMA designs, 
for ease of explanation, we have considered to focus on NPUs only. Then, as a case study, we consider 
wireless microphones as an example of NPU. The parameters of the OS-OFDMA are the same as in 
Section IIV-BI Before presenting the performance results we need to estimate the most realistic values of 
wireless microphones activity descriptors, i.e. average arrival rate and channel occupancy time. 

In the case of average NPU channel occupancy time we set it to a value between one and four hours, 
believing this represents a common activity time. More discussion is needed, however, on the arrival rate 
of the wireless microphones. Since the potential number of wireless microphones is location dependent, 
we have setup four different network scenarios, representing different places in the USA, see Table JIJ 
that differ in population density p and activity time 1/ p n . We assume that OS-OFDMA BS covers a 
fraction of the area of diameter L = 2 mi (for all scenarios) of the considered location, while the wireless 
microphones move in and out of the BS circular coverage with a certain speed v = 1.5 mph (for all 
scenarios). Then using a fluid flow model approximation ||63l we compute the average crossing rate 
of the wireless microphones to that area and translate it directly to an average arrival rate of wireless 
microphone on any of the subchannels. That is A n = phirLv, where \/h denotes number of inhabitants per 
one active wireless microphone in the considered location. Since the value of h is not known reliabljo we 
assume that one wireless microphone is present per 300 inhabitants (for all scenarios) and such wireless 
microphone is active for 10% of the time. Finally, U nmgx = m.a-K([-ir(L/2) 2 ph\, 1) in this case. We set 

6 The only credible report we were able to find was 1641 . The estimation using data present in this report was based on a 
simple calculation. According to |64. Sec. A2] there were 1924431 wireless microphone shipments in the European Union (EU) 
between 2002 and 2006, which translates to ~ 1 wireless microphone per 1000 EU inhabitants (assuming a constant level of 
wireless microphone shipments per year). Note that the value of 35,000-70,000 licensed wireless microphone operations in USA 
presented in (65] was not substantiated with any reference. 
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TABLE II 

Scenarios for the Analysis of NPU Activity Impact on the Performance of OS-OFDMA Designs 



Nickname 




"heavy urban" 


"urban" 




"light urban" 




"event" 


Location in US 




Los Angeles, 


CA 


Santa Barbara, 


CA 


Madison, WI 


Stapl 


;s Center, CA 


Users/mi 2 , p 




7452.7 




4708.2 




2701.0 




7452.7 


Activity time, l//i„ 


h 


1 




1 




1 




4 


^n,max 




8 




5 




3 
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the inter-arrival and departure times for other users as the same value as Section ITV-B I Also for analysis, 
we normalize all time parameters in the unit of five minutes as in Section IIV-BI 

The results are grouped separately for total average network throughput, CBR throughput and VBR 
throughput, see Fig. [3] Fig. @] and Fig. [51 respectively. We have chosen to vary number of CBR connections 
in all figures as a parameter, since in our model CBR connection is the most QoS sensitive and capacity 
demanding SU traffic class. First we immediately observe that SoNiBi implementation obtains the highest 
throughput. The larger the activity of the wireless microphones, the bigger the difference between SoNiBi 
and the remaining implementations - compare for example Fig. |3(c)| and Fig. |3(d)| Option S0N1B1 obtains 
a throughput around 3 Mbps, even in the "event" scenario. This is due to maximum utilization of the 
remaining channel capacity by subchannel notching and channel bonding. As the NPU activity increases 
all implementations reach almost zero throughput, while S0N1B1 still obtains reasonable performance. The 
worst performance is obtained for S0N0B0, while the S1N0B0 and S0N0B1 are in between the extreme 
cases. As the activity of the wireless microphones decreases all implementations start to converge in 
throughput - compare for example Fig. |4(a)| with Fig. |4(c)[ however the individual relation between 
the implementations stays the same. Then, we observe that all implementations except S0N1B1 obtain a 
similar throughput, irrespective of the network scenario. This proves that subcarrier notching promises 
to deliver most of the available capacity in the licensed bands. 

The worst situation, in terms of network scenario, is the "event" scenario. Due to long channel dwell 
time by NPU, i.e. four hours, the throughput for all implementations except S0N1B1 reaches zero. We also 
conclude that in scenarios where the activity of the wireless microphones is low, like in the "light urban" 



scenario, the users of systems based on OS-OFDMA are promised to obtain high QoS, see Fig. 3(c) 



Fig. |4(c)| and Fig. |5(c)[ irrespective of the implementation. 

A separate comment is needed for simulation verification of the results. In all cases S0N1B1 imple- 
mentation matches simulations perfectly, irrespective of the parameters selected. However, due to the 
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Fig. 3. Total throughput of different OS-OFDMA implementations for all network cases; Parameter setup and the ordering of 
the figures is the same as in Fig. [2] (a) "heavy urban", (b) "urban", (c) "light urban", (d) "event". 



approximations assumed for the remaining OS-OFDMA implementations, see again (|36l ), d40l ), the slight 
discrepancy is particularly visible for the scenarios with high NPU activity rate, compare for example 
Fig. |4(a)| with Fig. |4(c)| The results prove, on the other hand, that the developed model works very well 
for low PU activities, which is the typical case in real life PU occupancy statistics 0, |48l . Still, for 
each case study the relation between each OS-OFDMA implementation is well captured for any value 
of the parameters considered, while for the majority of the cases the mismatch between simulations and 
analysis is less than 10%. 
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Fig. 4. CBR throughput of different OS-OFDMA implementations for all network cases. Parameter setup and the ordering of 
the figures is the same as in Fig. [2] 



D. Impact of Two-stage Spectrum Sensing Options on OS-OFDMA Design 

The final experiment considers effect of sensing design parameters on the performance of OS-OFDMA 
designs. For this investigation we change the sensing time of coarse sensing, which has a direct effect 
on false alarm probability and makes a significant impact on the frequency of switching to fine sensing 
(and as a result on the throughput of the system). We keep other parameters the same as in the first 
experiment described in Section ITV-B1 except £/ njmax = 2, to see the effect of coarse sensing clearer, and 
U c ,ma,x = 10. We do not alter parameters of fine sensing phase, since we want to explore the benefit of 
two-stage sensing and the coarse sensing phase is a common element of every sensing method, including 
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Fig. 5. VBR throughput of different OS-OFDMA implementations for all network cases. Parameter setup and the ordering of 
the figures is the same as in Fig. [2] 



single stage sensing. 

Considering Rayleigh Channel with Additive White Gaussian Noise, we compute false alarm proba- 
bility, pio, and detection probability, pn, for an individual SU user as IJ351 Eq. (3)] and [35 . Eq. (4)], 
respectively. Then, according to IT351 Sec. III-B], we can derive p\\ as a function of p\o and r a for given 
average PU SNR and subchannel bandwidth b. Assuming collaborative sensing of all CBR and VBR 
users and OR logic for combining scheme, we compute system false alarm probability as 

<p a = l-(l-p 10 )^— +u —\ (23) 



and system detection probability, 5 a , as (|23l replacing pxo with p\\. In this evaluation we keep the 
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TABLE III 

Sensing Time and the Respective False Alarm Probability for the Experiment Introduced in Section HV-DI 
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Fig. 6. Throughput of OS-OFDMA implementations, (a) total, (b) CBR, and (c) VBR, as a function of coarse sensing time. 
Parameter setup is the same as in Fig. [2] except for U„ max = 2, (7c, max = 10 while r a and cf> a are given in Table [TTT] 



detection probability 5 a = 0.99 and change the sensing time for the coarse sensing such that r a = 
(0,4)ms. Note that r a = ms represents single stage sensing. Table ITTT1 presents calculated false alarm 
probability, 4> a , based on the assumed sensing time r a . 

The results are presented in Fig. [6] First interesting observation is that two stage sensing does not 
always provide better performance than single stage sensing. Comparing the total throughput at r a = 
with r a ~ (0.5, 1) ms of all OS-OFDMA options, two stage sensing shows worse performance than single 
stage sensing. This result is due to high probability of false alarm for this range of r a , see Table In] 
In addition, for this network setup, the throughput is maximized at r a = 2 ms and is larger than for 
T a = which confirms that two-stage sensing can indeed benefit all OS-OFDMA operations. This also 
confirms that the design choice of the IEEE 802.22 for spectrum sensing method was correct. When 
T a increases beyond the point for which the throughput is largest, the throughput of all OS-OFDMA 
implementations starts to rapidly decrease. This is due to the fact that the sensing overhead starts to 
dominate over potential improvement from decreased false alarm rate. This so called sensing-throughput 
tradeoff is in agreement with a similar investigation in the context of OS A ad hoc networks |[35l . Note 
that the relation between sensing time and the obtained throughput are the same when looking at total, 
CBR and VBR throughput, see the shapes of all curves in Fig. |6(a)[ Fig. |6(b)| and Fig. |6(c)| Also, the 
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different order in obtained throughput for each OS-OFDMA implementation are due to the specific OS- 
OFDMA options, not due to spectrum sensing parameters selected. For that, compare, e.g., the position 
of SiNqB in Fig. 6(a) with the position of the same implementation in Fig. |6(b)| For more on this aspect 



please refer to Section IIV-CI Finally, note that the simulation results match well with the analysis for all 
OS-OFDMA implementations. 

V. Conclusions 

In this paper we have proposed an analytical model that allows the comparison of different designs 
of OS-OFDMA. We have considered design options that include channel bonding, subchannel notching 
and two stage spectrum sensing. As a performance metric we have derived average throughput obtained 
by the secondary users of the spectrum. In the analysis we have included the inter-relations of different 
connection classes and priorities, like constant and variable bit rate traffic of the secondary users, and 
wideband and narrowband primary users. We concluded that OS-OFDMA design that allows the flexible 
bonding of channels and notching of subchannels currently occupied by the primary users, obtains the 
highest throughput in comparison to the designs that do not consider those options. As one of our 
numerical results show, the improvement reaches couple of hundred percent when the activity level of 
the different primary user types is very high. Also, as our investigation show, the two-stage spectrum 
sensing technique used in OS-OFDMA proves to increase the average network throughput, provided the 
probability of false alarm in the coarse sensing stage is low. 

Appendix A 
Derivation of G x (k\U x , X x ) and T x (j\U x , /j, x ) 

If the inter-arrival time has a negative exponential distribution with average arrival rate X x , the number 
of connections generated in a frame of length tt has a Poisson distribution. However, because we limit 
the maximum number of users to U x>mSjX , the number of users including newly generated connections, U x , 
cannot exceed U Xttnax . Considering this, we derive the probability of k new connections being generated 
in a frame, G x (k\U x ,X x ), as 



G x (k\U x , \ x ) 



(\ x t f ) k e x f , > „ jj jj 

g (Kt f ye-^f ^ k >0 }Ux = Uxniax . (24) 

i=k 



0, k < 0. 
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Note that the subscript x = {w, n, c, v} indicates the class of users, i.e., w for WPU, n for NPU, c for 
CBR, and v for VBR. 

Further, denoting the departure rate of each connection as /j, x , the probability that j connections are 
released from U x active connections during a frame of length tf, T x (j\U x , fj, x ,tf), can be calculated 
recursively as 

ft, 



T x (j\U x ,fx x ,tf) 



(25) 



V Ux^e-v^T^j -1\U X - l,tJt x> tf - t)dt, j > 0; 
e-u^tf, j = 0, 

which after some manipulation reduces to 

T x (j\U x ,v x ,t f ) = (U*\ e -u^t f(e ^t f _ iy (26) 

Since we only consider the connection release probability for users within the duration of a frame we 
abbreviate T x (j\U x , /i x ) = T x (j\U x , /J, x ,tf). 

Appendix B 
Derivation of Pri 3 (-) for S NiBi 

In Pri3(-) in (fTTT >. M a subchannels include M m subchannels that are occupied by PUs but mis-detected 
and Mo subchannels correctly detected as idle. Therefore, M a = M m + Mo, and thus to compute Pri3(-) 
we define a supporting probability, Pri4(M m ,Mo|C/ to , U n , S), which is the probability that the number 
of subchannels detected as idle correctly and falsely are Mo and M m , respectively. 

In design option So, a SU performs the coarse sensing for all subchannels first, and then, if the SU 
detects a PU on a subchannel, it immediately switches to fine sensing and senses all subchannels again 
with high sensing accuracy. Thus, for a condition 5 = (only coarse sensing is performed), the number 
of subchannels detected as idle must be the same as the number of all subchannels, i.e. M m + Mq = M. 
In other words, the case M m + Mo < M is impossible for the condition 5 = 0. For the condition 5 = 2 
(no idle subchannel is detected after all stages of sensing) only the case M m + Mo = is possible. For 
5 = 1, using the detection probability St and the false alarm probability 4>f in the fine sensing stage, 
we can derive the probability that M m busy subchannels are mis-detected and Mo idle subchannels are 
correctly detected for given U w WPUs and U n NPUs. Thus, 



Pru(- 



F 5 = 1; 

1, 5 = 0, M m + M = M, or 5 = 2, M m + M = 0; ( 27 ) 

0, 5 = 0, M m + M < M, or 5 = 2, M m + M > 0, 
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where F = (££) (1 - Sf)^'~ M " ( M ^) (1 - 0,)^"^-^ and M p = min(M, ^ + CW») 
is the number of subchannels actually occupied by PUs. Therefore, Pri3(-) can be derived as 

M a 

Pri 3 (-) = 5>r 14 (x,M a -x|[/ w ,[/ n ,S). (28) 

x=0 

Appendix C 
Derivation of Pri 5 (-) for S NiBi 

To derive Pris(-) we need to consider the following three cases. First, if the OS A network mis-detects 
existing PUs and correctly detects all idle subchannels in the coarse sensing stage, 5 = because OSA 
network will not advance to fine sensing. Second, if the OSA network detects at least one PU correctly 
or falsely in the coarse sensing stage and detects all subchannels as busy also correctly or falsely in the 
fine sensing stage, S = 2 because after coarse and fine sensing no idle subchannel is detected. Otherwise 
5 = 1. Thus, 

(l-S a ) M »(l-cf> a ) M - M », S = 0; 

,,i ) s (1 - (1 - 5 a ) M ?(l - <l>a) M - M ') Sf^f-^, 5 = 2; (29) 

1 - Pr 15 (5 = 0\U W , U n ) - Pr 15 (5 = 2\U W , U n ), 5=1. 



Pi- 



Appendix D 
Derivation of Pri 3 (-) for S N Bi 

For the case SoNoBi, because the spectrum is utilized based on the unit of bandwidth of one channel 
rather than one subchannel, we introduce the probability of the number of channels detected as idle, X a , to 
compute the probability of the number of subchannels detected as idle, M a . Denoting Pr2o {X a \ U w ,U n ,S) 
as the probability that the number of channels detected as idle is X a given U w WPUs and U n NPUs in 
the 5 sensing stage, (|28T ) is modified as 

a |Pr 20 (X o |C/ w ,C/ n ,5), M a = X a Y; 
PriaCO = \ (30) 

I 0, otherwise. 

There are two required conditions for Pr2o( - ) hi (|3Qb . First, the sum of the channels detected as idle and 
the number of WPUs cannot be greater than the total number of channels, i.e. X a + U w < X. Second, 
the number of NPUs cannot be greater than the total number of subchannels that are not occupied by the 
WPUs, i.e. U n < (X — U w )\Y/l n \. We denote l c as an indicator of those conditions, defining l c = 1 
for X a + U t < X or U n < (X - U w ) \Y/l n \ and l c = 0, otherwise. 
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Next, we compute Pr2o(-) under the feasible conditions considering different sensing stages. Under 
the condition that no channel is detected as idle after the fine sensing, i.e. 5 = 2, the number of channels 
detected as idle is X a = under the assumption of perfect PU detection. On the other hand, when all 
channels are detected as idle in the coarse sensing, i.e. 5 = 0, only X a = X is possible. The analysis for 
the condition that the SU performs fine sensing and detects idle channels, i.e. 5 = 1, is not easy to derive 
directly because the number of channels detected as idle depends on the position of NPUs in the spectrum 
as well as the number of NPUs. Thus, defining the number of channels actually occupied by the NPU as 
X n , we deconstruct Pr2o(-) for 5 = 1 into two components: the conditional probability of X n for given 
number of PUs and the sensing stage, denoted as ~Pi2i(X n \U w , U n , 5), and the conditional probability 
of the number of channels detected as idle for a given X n , denoted as P^^al^n, U w , U n , 5). Then, 

Y, x ?T2i{X n \U w , U n , S)Pv 22 (X a \X n , U w , U n , 5), l c = 1, S = 2; 

l c = land 5 = 0, 
rv ; > n () = <|l, X a = XorS = 2, (31) 

X a = 0; 
0, otherwise, 

where X = {X n \\(U n l n /Y~\ < X n < mm(U n ,X — Ut — X a )} because X n is the smallest when all 
NPUs are located on adjacent subchannels, i.e. \(U n l n )/Y~\, and the largest when all NPUs are located 
in different channels separated as far as possible, i.e. mm(U n ,X — Ut — X a ). 

To compute Pr2i(-) in d3"Tl ), we assume that any NPU can appear on any subchannel with equal 
probability, and that false alarm can occur uniformly over all idle subchannels. Then, we introduce a 
supporting function f s (k,x,r) denoting the number of possibilities that k items are distributed over 
exactly x bins each of which has a capacity of r items. 

To derive f s (k,x,r) first we introduce the supporting variable, ij - the number of items in j-th bin 
where j € {1, • • • , x}. Then, there can be (H possible distributions for the j-th bin, and thus for x bins 
there can be Y?j=i CO ' ' ' Y^j=i O possibilities. Because there should be no bin empty, ij > 1. Also, 
each of the bins has a capacity of r items and the total number of items cannot be greater than k, and 
therefore ij < min(r, A;). In addition, if the number of bins is less than the number of items or equal to 
zero, there is no way to fill all x bins. Thus 

f s (k,x,r) = \ El(k < x < r) E}U ^ ' ' ' E5Ll G) ' ° <X< k] (32) 

I 0, otherwise, 
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where I(k,x,r) = {h,--- ,ix\YTj=\ i j = k,h,- ■ ■ ,i x e {I,- ■ ■ ,min(r, k)}}. 

Then Pr2i(-) can be computed by dividing the number of possible events that U n NPUs are located 
on exactly X n channels each of which can have at maximum L^AnJ NPUs, i.e. f s (U n ,X n , \Y/l n \), by 
the number of all possible events that U n NPUs appear on X — U w channels, i.e. (' ~ j^H '"■'). Note 
that if there is no NPU, i.e. U n = 0, then X n should be zero. Considering the possible case of selecting 
X n NPU channels from a total of X — U w channels, i.e. ( x "), we derive Pr2i(-) as 

1, x n = o,u n = o, 

Pr 2 i(-) = < { x-v*. )MUntXnt[Y/lni) . (33) 

^ — ", {x - Uwn Y/i n i\ , otherwise. 

Now we present the derivation of Pr22(-) i n *EB- If there exist U w and X n channels that are occupied 
by WPUs and NPUs, respectively, and X a channels are correctly detected as idle, the remaining X — 
U w — X n — X a channels must be falsely detected as busy. Considering the number of events of selecting 
X a channels from X — U w — X n idle channels, we derive Pr22</) as 

Pr 22 (-) = (* " \~ Xn ) (1 - (1 - ^rf'^~ X ^(l - 0/)^. (34) 

Appendix E 
Derivation of Pr 23 (-) for for S N B! 

Because for SoNoBi the spectrum is utilized based on the unit of bandwidth of one channel, denoting 
a supporting probability Pr 24 (xjP \XJf~V , U$ , U^~ l \ U^ , C^* _1) , S®, S^A as the conditional 
probability of the number of channels, we derive Pr23(-) as follows 

L lPr M (X^\xt 1 \uS\ut 1 \v^\ut 1 \S^S^-^), M a = X a Y; 
Pr 23 (0 = I (35) 

I 0, otherwise. 

To reduce the complexity of calculating Pr24(-) we use an approximation that if the number of PUs 
is the same at times t and t — 1, the positions of NPUs at times t and t — 1 are also the same. This 
approximation is valid when the PU activity is not high. With this approximation we ignore the case that 
a certain number of NPUs disappear while at the same time the same number of new NPUs appear, but 
at different locations. Also we assume that if the number of WPUs and NPUs have changed, the number 
of channels detected as idle is independent of the number of WPUs and NPUs at time t — 1, which means 
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Pr 



2H 



that Pr24(-) 2 Pr2o(-)- Considering these approximations we derive Pr24(-) as 

1 Tj{t) _ TT(t-l) Tj(t) _ rr(t-l) X (t) _ v{t-l). 

o, t/w = ^-D, u® = ut'lxp + jr<*-i) ; (36) 

Pf20 (-^a I C^u j f^n > <5) > otherwise, 
where Pr2o( - ) is defined in Appendix ID1 as (DTK 

Appendix F 
Derivation of Pri 5 (-) for S N Bi 

To modify Pri5(-) in d29l ), under the approximation of the perfect detection, depending on the number 
of subchannels occupied by PUs, i.e. M p , we consider three cases. In the first case, M p = M, all 
subchannels are occupied by PUs. Then, because of the perfect detection approximation S = 2 is the 
only feasible condition. In the second case, M p = 0, no PU appears on the subchannels. In this case, the 
false alarm probability of each sensing stage affects the probability Pris(-). In the last case, < M p < M, 
the probability of S = is zero and the false alarm probability of the coarse sensing does not affect the 
performance. Thus 
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0< M p <M,S = 1; 
< M p < M, S = 2; 



0, 



otherwise. 



Appendix G 
Derivation of Pri 5 (-) for SiN B 

Following the same convention as for the derivation of © in Section IlII-AI denoting Pt2g(S, U w ,U n ) 
as the joint probability of the number of PUs and the sensing stage, we derive Pris(-) as 

p M a Pr 2 e(S,U w ,U n ) 

PTl7{U w ,U n ) 

Because expression Pr26(-) is the steady state probability, denoting the state transition probability as 
Pr 27 (S®, U$, 1$ \S { - t - 1 \ t4f _1) , U% _1)N ) , Pr 2 e(-) can be obtained by solving the Markov chain with 
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a state {S, U w , U n }, such that Pr 2 6(-) = Ssc-d P r 26(-) Pr 27(') and ^5 Pr 26(') = 1. where S is the set 
of all possible S, U w , and U n , and 5^ -1 ^ is the set of the same variables at time t — 1. Then Pr27(-) is 
computed as 

PrarC-) ± Pr 25 (&\u£\U<!\S^ 1 \ut 1) ,Ut 1) ) Prie (ltf\ t^l^, itf" 1 ') • (39) 

To compute Pri2(-) in d22l) and Pris(-) in (l38l) . we need to derive Pr2s(-)- For the derivation of Pr2s(-)> 
to avoid prohibitive complexity of calculation, we consider two specific cases. The first case is that there 
is no change for the number of WPUs and NPUs between time t — 1 and t. This case is expected to 
happen most frequently because the PU generally has longer inter-arrival time and packet length than 
the frame size of the SUs. The second case is that the numbers of WPUs and NPUs are reduced as time 
goes from t — 1 to t. If the PU activity is not too high, this case happens when some of the existing PUs 
disappear and no new PU enters into the subchannels. Then subchannels detected as idle at time t — 1 
still remains idle even at time t, which can affect the sensing stage significantly. For the other cases, we 
use a general approach. Depending on the changes of the numbers of the PUs, Pr2s(-) is derived as 

'Pr 25a p), U$ =ut X \u® =ut X \ 
Pr 25 (-) « <! Pr 25b (D) , U® < Ut l \ U& < Ut 1] , (40) 

Pf25c (D) , otherwise, 

where D = {S®\U$, U^ , S^'^, U^ , U^ 1] }. To analyze the first case in (HI, where there is no 
change in the number of WPUs and NPUs between time t — 1 and t, we reuse the assumption that 
there is no change in the positions at subchannels of the NPUs for the same number of NPUs at times 
t and t—1, see Section ITlI-Bi With this approximation, under the condition that there exists at least one 
channel at time t — 1, i.e. S^' -1 ) < 2, there should also exist at least one idle channel at time t, i.e. 
gw < 2. Thus, if no false alarm occurs, only coarse sensing is performed, i.e. £W = and otherwise 
S^' = 1. Next, for the condition that no channel is detected as idle at time t — 1, i.e. S^ 1 ' = 2, there 
can be a large number of possible events depending on the positions of NPUs and false alarms. Thus, 
instead of considering all the possible events, we consider two cases and apply approximations for each 
one. The first case is that there exists a small number of NPUs at time t — 1 so that the NPUs can not 
occupy all X — Uw channels, i.e. f s (Un ,X — U w , |_^V^nJ ) = 0, where / s (-) is defined in 
Appendix [D] as (|32~1) . In this case, because there exists at least one idle channel, we approximate that 
the SU performs fine sensing and detects the idle channel at time t, i.e. S® = 1. In contrast, for the 
second case where there exists enough NPUs such that they are spread over all X — U w ~ channels, 
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i.e. f s (U„. ,X — Uw , [Y/l n \ J > 0, we approximate that all channels are detected as busy at time 
t, i.e. S® = 2. Considering all these sub-cases, we derive Pr25 a (-) as 



Pr 2 5a(-) 



(1-4 
i-(i 

i, 

o, 
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s (t) = 0> s (t-i) < 2 . 

5 (t) = 2> j(*-i) > o, or S® = 1, J^- 1 ) 
otherwise, 



(41) 



where J^ = / 8 [U^',X - L4 , iTAnJ )• Next, we analyze the second case in dH, i.e. U$ < U, 
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r(*-l) 
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and Un < Un , which implies that the number of PUs decreases. Similar to the first case in (l40l ). 
if there exist idle channels at time t — 1, i.e. S^* -1 ) < 2, there must also exist idle channels at time t, 
i.e. S® < 2. If a false alarm does not occur on any of the Y subchannels in a channel utilized for data 
transmission, the SU performs only coarse sensing, i.e. S^' = 0, and otherwise, S^' = 1. On the other 
hand, for the condition that there is no idle channel at time t — 1, i.e. S^*" 1 ) = 2, the probability of the 
sensing stage at time t is calculated by dividing the number of the cases that U n NPUs occupy exactly 



all X - U$ channels, f s (u$ , X 



U%\ [Y/l r 



, by the number of all possible cases 



f (x-uW)lY/l n ] 



)■ 
Note that when S^ -1 ) = 2, it is not possible to perform only the coarse sensing at time t, because the 

channel that is going to be utilized for data transmission is not determined. In other words, the probability 

of S^ 4-1 ) = under 5 (f_1) = 2 equals zero. As a result, the probability of S^"^ = 1 under S^' 1 ^ = 2 

equals to one minus the probability of S^ -1 ) = 2 under S^' -1 ) = 2. Thus, we derive Pr 2 5b(-) as 



Pr 2 5fe(-) 



(l"0a) y , 
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o, 
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f(X-U^')yY/l„\-y 



5W=0,5('- 1 )<2; 
5(*) = l ) 5(*-i)<2; 
5(t) =2) 5(*-i) = 2; 

5 (i) = 1)5 (t-i) = 2; 
otherwise. 



(42) 



Finally, we compute Pr25 C (-) in d4Q~b . In this case, again to simplify the computation, we ignore the 
effect of the sensing stage at time t — 1 and only focus on the sensing stage at time t. Then, we compute the 
probability that the SU performs only coarse sensing, i.e. S^ = 0, by dividing the number of cases that 



at least one channel remains idle and no false alarm occurs on that channel, (1 — <fi a ) 



by all possible cases of all channels remaining idle, 



f(X-TJW)\Y/l n 

^ uk t] 



Y((X-UW-l)[Y/l„} 
\ UP 



)• 



On the other hand, the probability 
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that the SU detects no idle channel, i.e. S® = 2, is calculated by considering the case that U n NPUs 
occupy exactly all X — U w channels, i.e. f s (U^,X — U w , \Y/l n \). Note that for the case that there 
appears too many NPUs on the spectrum so that there is no possibility that even one channel cannot 
remain idle, only S = 2 is feasible. This case can be denoted as V = (X — U w — l)\Y/l n \, because 
X — U w channels are available for the NPU and each channel can have a maximum L^AnJ NPUs. 



Considering all these cases, we have 



Pr 25c (-) 



<t a. \Y(( x - u ™ ) - i nY/i„i\ 

f(X-U^)[y/lni\ 



S(*) = 0, UX' < V; 



J^l alt) - o rrW 



((X-u£ ] )YY/l n \\ ' 








U — A, u n ^ V 


1 - Pr 25c {S® -- 


= o|-)- 


- Pr 25c (SW = 


= 2|"). 


S(*) = l, ujp < V 


1, 








S(*) = 2, UJp > V 


0, 








otherwise. 



(43) 
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